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Electronic structure calculations combining the local-density approximation with an exact diago- 
nalization of the Anderson impurity model show an intermediate 5/^-5/^-valence ground state and 
delocalization of the 5/^ multiplet of the Pu atom 5/-shell in PuCoIns, PuCoGas, and (5-Pu. The 
5/-local magnetic moment is compensated by a moment formed in the surrounding cloud of conduc- 
tion electrons. For PuCoGas and 5-Pu the compensation is complete and the Anderson impurity 
ground state is a singlet. For PuCoIns the compensation is partial and the Pu ground state is 
magnetic. We suggest that the unconventional d-wave superconductivity is likely mediated by the 
5/-states antiferromagnetic fluctuations in PuCoIns, and by valence fluctuations in PuCoGag. 

PACS numbers: 74.70.Tx, 74.45. +c, 74.20.Mn, 74.20.Pq 



Providing a consistent description of correlation ef- 
fects in the electronic structure of elemental actinides 
and their compounds is a complex problem due to 
the interplay between the localized and the itinerant 
nature of the 5f electrons. It is commonly accepted 
that 5/-electrons in light actinides form rather broad 
conduction bands whereas for the heavy actinides the 
5/ states are atomic-like. Johansson [l| described this 
situation as a "Mott transition in the 5/-electron sub- 
system" taking place between Pu and Am when mov- 
ing along the Periodic Table. Katsnelson et al. Q 
linked the broadening of the 5/ band to the "atomic 
collapse" characterizing the transformation from the 
high-temperature expanded and the low-temperature 
compressed phases of Pu. 

A quantitative description of the Mott transition 
in actinides Q was obtained by the dynamical mean- 
field theory (DMFT) [3| more than 20 years after the 
concept was formulated. Further DMFT studies sug- 
gested an intermediate- valence nature of the Pu-atom 
5/ shell and provided justification for the experi- 
mentally proved absence of magnetism in (5-Pu [g] . 

The intermediate-valence and nonmagnetic char- 
acter of the 5/ shell can play an important role 
in stabilizing the superconducting state exhibited by 
PuCoGa5 below a critical temperature of 18.5 
The unconventional character of supercon- 
ductivity in this compound is now generally accepted 
but the microscopic mechanism responsible for elec- 
tron pairing remains unknown. The d-wave symme- 
try of the superconducting gap in PuCoGas has been 
proven by point-contact spectroscopy experiments [l^l 
that also provided the first spectroscopic measure- 



ments of the gap amplitude and its temperature de- 
pendence. 

Recently, superconductivity has been discovered 
also in PuCoIns [13, with = 2.5 K. The exper- 



imental studies of this compound were immediately 
followed by conventional density functional theory 
(DFT) calculations in the local-density generalized- 
gradient approximation (LDA/GGA) [12|, |l3[. Keep- 
ing in mind a well known failure of DFT in the case 
of (5-Pu 0, it can be expected that LDA/GGA docs 
not provide an accurate description of the electronic 
structure for this strongly correlated material. A few 
static mean-field correlated band theory calculations 
were also performed [13, [3 , making use of different 
flavors of the LDA/GGA plus Coulomb's U (LBA+U) 
method. While being an improvement over the con- 
ventional band theory, the LDA(GGA)-|-J7 falls short 
in describing the itinerant-to- localized crossover of the 
5/ manifold in S-Fu [3 and PuCoCag 

Here, we report electronic structure calculations of 
PuCoIns, PuCoGas and 6-Pu performed by combin- 
ing LDA with the exact diagonalization (ED) [l5| of 
a discretized single- impurity Anderson model [l6| . In 
this approach, the band structure obtained by the rel- 
ativistic version of the full-potential linearized aug- 
mented plane wave method (FP-LAPW) [T^ IS con- 
sistently extended to account for the full structure of 
the 5/-orbital atomic multiplets and their hybridiza- 
tion with the conduction bands [l^ ■ 

The starting point of our a£ 
band Hubbard Hamiltonian [li 



jroach is the multi- 
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one-particle Hamiltonian found from ab initio elec- 
tronic structure calculations of a periodic crystal; iJ™* 
is the on-site Coulomb interaction [l^ describing the 
/-electron correlation. We assume that electron inter- 
actions in the s, p, and d shells are well approximated 
in DFT. 

The effects of the interaction Hamiltonian i/'"' 
on the electronic structure are described by a k- 
independent one- particle selfenergy S(z), where z is a 
(complex) energy. The selfenergy is constructed with 
the aid of an auxiliary impurity model describing the 
complete seven-orbital 5/ shell. This multi-orbital im- 
purity model includes the full spherically symmetric 
Coulomb interaction, the spin-orbit coupling (SOC), 
and the crystal field (CF). The corresponding Hamil- 
tonian can be written as [l6l | 
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where fXia creates an electron in the 5/ shell and 6j,jg. 
creates an electron in the "bath" that consists of those 
host-band states that hybridize with the impurity 5/ 
shell. The energy position e/ of the impurity level, and 
the bath energies are measured from the chemical 
potential /i. The parameter ^ specifies the strength 
of the SOC and Acf is the crystal-field potential at 
the impurity. The parameter matrices describe 
the hybridization between the 5/ states and the bath 
orbitals at energy e'^ . 

The band Lanczos method is employed to find 
the lowest-lying eigenstates of the many-body Hamil- 
tonian -ffimp and to calculate the one-particle Green's 
function [Gimp{z)\'^^, in the subspace of the / or- 
bitals at low temperature (fesT = 1/500 eV). The 
self-energy is then obtained from the in- 

verse of the Green's- function matrix [Gimp]- 

Once the self-energy is known, the local Green's 
function G{z) for the electrons in the solid, 



[G{z)] 



71 72 



^|^^d3fc[z+M-if,DA(k)-E(z)];;^, 

(2) 

is calculated in a single-site approximation as given 
in Then, with the aid of the local Green's 

function G(z), we evaluate the occupation matrix 

'7172 — TT 71 72- The matrix 7^7l72 



is used to construct an effective LDA-|-t/ potential Vu, 
which is inserted into Kohn-Sham-like equations: 

[-V^ + FldaW + Vu+ C(1 ■ s)]<(r) = e^<(r). 

(3) 

These equations are iteratively solved until self- 
consistency over the charge density is reached. In each 
iteration, a new Green's function Glda(-z) [which cor- 
responds to G{z) from Eq.(l2]) with the self-energy S 
set to zero] , and a new value of the 5/-shell occupa- 
tion are obtained from the solution of Eq. ([3|) . Subse- 
quently, a new self-energy E(z) corresponding to the 
updated 5/-shell occupation is constructed. Finally, 
the next iteration is started by evaluating the new 
local Green's function, Eq.®. 
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FIG. 1. (Color online) The Pu atom LDA j = 5/2, 7/2 
projected DOS, and LDA hybridization function A(e) = 
-^ImTr[G'^ [e + iS)]. The inset shows the PuCoIns crys- 
tal structure. 

In order to determine the bath parameters 
and , we assume that the LDA represents the 
non-interacting model. We then associate the LDA 
Green's function Glda(2:) with the Hamiltonian of 
Eq. ([T]) when the coefficients of the Coulomb interac- 
tion matrix are set to zero {Umm'm"m"' — 0). The hy- 
bridization function A(e) is then estimated as A(e) = 
-i ImTr[GLDA(e + i6)]. The curve obtained for A(e) 
is shown in Fig. [1] together with the j = 5/2,7/2- 
projected LDA densities of the /-states. The results 
also show that the hybridization matrix is, to a good 
approximation, diagonal in the {j,jz} representation. 
Thus, we assume the first and fourth terms in the im- 
purity model, Eq. ([T]), to be diagonal in {j,jz}, so that 
we only need to specify one bath state (six orbitals) 
with 6*^^5/2 s-^d V^*|!^/2' ^^"^ another bath state (eight 
orbitals) with £j^y2 a-^d V^^^y^. Assuming that the 
most important hybridization is the one occurring in 
the numerical values of the bath 



n. 



the vicinity of Ep, 

'^5/2,7/2 



parameters V^^/^l^^^ are found from the relation [20| 
-A{e)/Nf integrated over the 
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energy interval, Ep — 0.5 cV < e < Ep + 0.5 eV, with 
Nf = 6 for j = 5/2 and Nf = 8 for j = 7/2. The 
bath-state energies ji shown in Tabic |T] are ad- 
justed to approximately reproduce the LDA 5/-state 

7/2 



occupations n^"^ and 



TABLE I. 5/-states occupations nl"^ and w^^, and bath 
state parameters £5/27/2 (^V), (^V) for Pu-atom 

in PuCoIns, PuCoGas, and 5-Pu from LDA calculations. 



Material 


5/2 


7/2 


4^' 




7/2 
^1 




PuCoIns 


4.78 


0.39 


0.36 


0.21 


-0.25 


0.25 


PuCoGas 


4.38 


0.76 


0.25 


0.29 


-0.07 


0.34 


5-Pn 


4.16 


0.85 


0.33 


0.27 


-0.01 


0.36 



In the calculations we used an in-house implemen- 
tation Ell, [H of the 



, FP-LAPW method that in- 

cludes both scalar-rclativistic and spin-orbit coupling 
effects. The calculations were carried out assuming 
a paramagnetic state with crystal structure param- 
eters for PuCoIns, PuCoGas, and S-Pu taken from 
Rcfs. [llj, |23|, |2J], respectively. The Slater integrals 
were chosen as Fq = 4.0 cV, and F2 = 7.76 eV, F4 — 

They corresponds 
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5.05 eV, and Fg = 3.07 eV 
to commonly accepted values for Coulomb's U = 4.0 
eV and exchange J = 0.64 eV. The SOC param- 
eters ^ = 0.28 eV for PuCoIns and PuCoGas and 
0.29 eV for S-Pu were determined from LDA calcula- 
tions. CP effects were found to be negligible and Acf 
was set to zero. For the double-counting term enter- 
ing the definition of the LDA+U potential, Vu, we 
have adopted the fully-localized (or atomic-like) limit 
(FLL) Vdc = U{nf -1/2) -J{nf-l)/2. Furthermore, 
we set the radii of the atomic spheres to 3.1 a.u. (Pu), 
2.3 a.u. (Co), 2.3 a.u. (Ga), and 2.5 a.u. (In). The 
parameter Rp^ x i^max = 10.54 determined the basis 
set size, and the Brillouin zone (BZ) sampling was 
performed with 1152 k points. The self-consistent 
procedure defined by Eqs. (P)-© was repeated un- 
til the convergence of the 5/-manifold occupation nj 
was better than 0.01. 

Wc arc now ready to discuss the solution of Eq.lJT]). 
For PuCoIns, the ground state of the cluster formed 
by the 5/ shell and the bath is given by a superposi- 
tion of a magnetic sextet (23%) and a non-magnetic 
singlet (77%), with occupation numbers (n/) = 5.40 
in the / shell and {rn,ath) ~ 8.40 in the bath states. 
This ground state is not a singlet and carries a non- 
zero magnetic moment. For the 5/ shell alone, the ex- 
pectation values of the spin (5*/), orbital {Lf) and to- 
tal (J/) angular moments can be calculated as {X'j) ~ 
XfiXf + 1) {Xf = Sf,Lf,Jf), giving Sf = 2.27, 
Lf ^ 3.90, and ,// = 2.09. The individual compo- 
nents of the moments vanish, (Sj) = (Lj) = 0, unless 
the symmetry is broken by an external magnetic field. 



In the case of PuCoGas, on the other hand, the 
hybridized ground state of the impurity is a non- 
magnetic singlet with all angular moments of the 5/- 
bath cluster equal to zero {S = L = J = 0). It consists 
of {nf) = 5.30 / states and (ribath) = 8.70 bath states. 
In a pictorial way, we can imagine that the magnetic 
moment of the 5/ shell (for which we get Sf = 2.18, 
Lf = 4.05, Jf = 2.43) is completely compensated by 
the moment carried by the electrons in the conduction 
band. As the value of the 5/ magnetic moment fluctu- 
ates in time, because of the intermediate- valence elec- 
tronic configuration, this compensation must be un- 
derstood as dynamical in nature. The same situation 
is realized in S-Pu {Sf = 2.11, Lf = 4.21, Jf = 2.62), 
whose ground state is found to be a nonmagnetic sin- 
glet with (nf) = 5.21 and {math) = 8.79. 

The 5/-orbital density of states (DOS) obtained 
from Eq. ([2|) for the three investigated compounds 
is shown in Fig. [2l Below the Fermi energy Ep the 
DOS exhibits the three-peak structure typical for Pu 
and for a number of its compounds, and its shape is 
in good agreement with experimental photocmission 
spectra. It can be noticed that the multiplets for the 
atomic configuration (/^ — s- transition, lying 
closer to Ep) are better resolved than for the /"^ part 
of the spectrum (/^ — >• transition). 

Comparison with previous LDA-|-Hubbard-I (HIA) 
calculation for S-Pu [l3l, and PuCoGas 01 shows 
that the three-peak manifold lying above 2 eV bind- 
ing energy has a slight upright shift towards Ep. At 
binding energies around 4 eV, the LDA-I-HIA peaks 
are substantially modified, and in the LDA-fED cal- 
culations they are spread over a ~ 3 cV energy in- 
terval. These changes in the DOS are induced by the 
hybridization and suggest partial delocalization of the 
multiplet. This is a situation suggested first by 
Hanzawa [27l | in intermediate- valence rare-earth com- 
pounds such as SmS or SmBg, where fluctuations oc- 
cur between two atomic-like 4/ conflgurations. Here, 
the 5/ states remain localized for the configuration 
but become itinerant for the one. 

As the many-body resonances lying closer to the 
Fermi energy are produced by multiplet tran- 

sitions, they are in a way analogs to the Racah peaks, 
specific transitions between Racah multiplets [28| of 
/" and /"^^. On the other hand, these structures 
determine the metallic character of the investigated 
materials that can therefore be considered as a realiza- 
tion of a Racah metal, situated between the two limit- 
ing cases represented by fully localized intermediate- 
valence rare-earth compounds and metallic systems 
(e.g., nickel) with a non- integer number of d electrons. 

Both PuCoGas and S-Pu display a temperature- 
independent magnetic susceptibility at low tempera- 
tures [2^ . Analogous to the intermediate- valence 
rare-earth compounds f30j . the magnetic susceptibil- 
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FIG. 2. (Color online) /-electron density of states (DOS, 
j = 5/2, 7/2 projected) for the Pu atom in PuCoIns (a), 
PuCoGas (b) and 5-Pu (c). 



at Ep, and 
{Tt[N{Ef){1 

Tfc = 72 mcV (- 850 K) for PuCoGag and T/e = 63 
mcV (^ 750 K) for 5-Pn. Since Tfc is high, x remains 
constant for T <C Tfc, as observed experimentally for 
PuCoGas and S-Pu. The situation is different in the 
case of PuCoIns where the ground state of the impu- 
rity is not a pure singlet due to weaker hybridization. 
Consequently, the temperature dependence of x is ex- 
pected to be more pronounced. 

The electronic specific-heat coefficient can be esti- 
mated as 7 = ^klTv[NiEF){l - ^)U=o]. For 
(5-Pu, we get w 44 mJ mol~^, in very good agree- 
ment with experimental data. For PuCoGas, we get 
« 43 mJ mol~^ which is smaller than the experi- 
mental value of 80-100 mJ mol~^. For PuCoIns, 
the estimated 7 value of ~ 52 mJ mol~^ is even 
further away from the experimental value of 180 mJ 
K^^ mol^^. In this case, it is difficult to obtain an ac- 
curate value for 7 due to the sharp DOS peak in the 
vicinity of Ep (see Fig. [5]). When taken right at the 
DOS peak position, the 7 value of 95 mJ mol~^ 
is obtained. Also, note that a possible enhancement of 

' 1.0 ' ' 
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ity is anticipated to behave as x 1/(2^ + ^^c); where 
the temperature Tfc describes fluctuations between 
the 5/ and conduction band electron states. Tfc cor- 
responds indeed to the broadening of the quasiparticle 
resonance near Ep due to valence fluctuations |31| . As 
the ground state of the impurity is a singlet, we esti- 
mate Tfc using a renormalized perturbation theory of 
the Anderson model Tfc = ~^Z[A{Ep)/Nf], 
where [A{Ep)/Nf] is the hybridization per orbital 



FIG. 3. (Color online) (Top) The band structure with /- 
weight fatbands for PuCoIns , and (bottom) the Fermi sur- 
face of PuCoGas and PuCoIus obtained from LDA-I-ED 
calculations. The shade of colors encodes the size of the 
energy gradient. 

Figure [3] shows the band structure and the eorre- 
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spending Fermi Surface (FS) for PuCoIns, calculated 
from the solutions of Eq. ([3|) , which represents an ex- 
tended LDA-I-J7 static-mean-field band structure with 
the 5/-states occupation matrix obtained from the lo- 
cal impurity Greens function Eq.Q. For comparison, 
Fig. pshows also the FS for PuCoGag (Fig. S2 of 
Ref. [10[). Close similarities in the band structure of 
the two compounds are immediately apparent. Both 
are compensated multiband metals, as the Fe-based 
superconductors, and for both materials the / bands 
move away from the Fermi level when the Coulomb- 
U is included, as can be seen by examining the /- 
weighted fatbands. The Fermi surfaces are composed 
by four sheets (1-4), one that is hole- like (FS-1) and 
three that are electron-like (FS-2,3,4). The Fermi ve- 
locities ratio {vl^y)^ / {vD^ of 1.54 for PuCoIns, and 
1.55 for PuCoGas are calculated in reasonable agree- 
ment with the experimental anisotropy ratio of the 
critical field iJc2, 2 — 2.3 for PuCoIns, and indicate a 
two-dimensional character of the electronic structure. 

DFT electronic structure calculations for Pu-based 
115 material have recently been reported by Ronning 
et al. 13 1 and Zhu et al. [l2j Their analysis of the 



DFT band structure and FS (see, e.g.. Figs. 3 and 4 of 
Ref. [lH) indicated two possible superconducting gap 
symmetries, the so-called s± and d^2_y2, which cor- 
respond to a pairing potential peaked at the (tt, tt, 0) 
reciprocal lattice position. The conclusion was drawn 
that for Pu-based "115" superconductors, the s± or- 
der parameter is more likely that the d^^-yi one. This 
is in contradiction with point-contact spectroscopy re- 
sults [l^ showing a zero-bias conductance anomaly 
that is not expected for s± gap symmetry 32 1. 



The presence of a 5/ local moment dynamically 
compensated by the surrounding conduction electrons 
together with the p-f^ intermediate-valence ground 
state in PuCoGas and PuCoIns opens various pos- 
sibilities for unconventional superconductivity. In 
PuCoIns the Pu /-shell local moment is not fully com- 
pensated and superconductivity could be related to an 
antiferromagnetic quantum critical point [ill [33j . On 
the other hand, in PuCoGas the ground state is a 
singlet and it seems more plausible that superconduc- 
tivity results from a valence instability, as in heavy- 
fermion superconductors (s^ l . 
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